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Abstract 



Of 

, Particle creation during inflation is considered. It could be important for 

species whose interaction is of gravitational strength or weaker. A complete 
but economical formalism is given for spin-zero and spin-half particles, and the 
particle abundance is estimated on the assumption that the particle mass in 
the early universe is of order the Hubble parameter H. It is roughly the same 
for both spins, and it is argued that the same estimate should hold for higher 
spin particles in particular the gravitino. The abundance is bigger than that 
from the usual particle collision mechanism if the inflationary energy scale is 
■ of order 10 16 GeV, but not if it is much lower. 

g> ■ 

I. INTRODUCTION 

a.: 

According to present ideas, particles species can be created in the very early universe by a 
number of mechanisms |D],|2|]. The most familiar, which applies to particles of any spin, is the 



collision or decay of other particles. A number of other reasonably familiar mechanisms come 
into play for the case of a spin zero particle, which corresponds to the oscillation of a scalar 
field. If the field is initially displaced from its vacuum value it will oscillate homogeneously, a 
familiar example being the oscillation of the inflaton field after inflation. The homogeneous 
oscillation can in turn cause further particle production, either through the decay of the 
initial particles as in the original discussion of reheating after inflation, or through collective 
effects as in the more modern discussions of 'preheating'. Alternatively, the scalar field can 
be radiated from topological defects, as in the example of axion production from axionic 
strings. 

The present paper explores a mechanism which is different from any of the above, namely 
the production of particles from the vacuum fluctuation during inflation. We estimate the 
abundance of particles arising from the vacuum fluctuation during inflation, and compare 
it with the abundance from the particle collision mechanism and with the cosmological 
constraints. 

In contrast with earlier work, we take seriously the fact that the effective mass m(t) in 
the early universe is not expected to be the same as the true mass m which is relevant at 
the present day, because of the effect of interactions. On the usual assumption that we are 
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dealing with a supersymmetric theory, it has long been recognised that even an interaction 
of gravitational strength will typically give a spin zero particle a mass of order H during 
inflation if H is bigger than the true mass More recently it has been emphasised || 

that the same result is to be expected even after inflation. Although it has not so far been 
stated explicitly in the literature, one may expect this result to apply equally to particles of 
any spin. 

This contribution of order H to the mass can be evaded if the mass is protected by a 
symmetry. For example gauge bosons cannot acquire mass unless the symmetry is sponta- 
neously broken, and Peccei-Quinn symmetry should keep the axion mass small, but neither 
of these cases is of interest in the present context (gauge bosons are too short-lived and the 
axion density will be dominated by other production mechanisms). 

During inflation there are more possibilities for evading the mass of order H . A number 
are available for spin zero, as reviewed for example in for the particular case of the inflaton 
field. Another case is that of the gravitino; its mass is proportional to \ W\ where W is the 
superpotential, and in some models of inflation |]5]|| it is much less than H during inflation. 
However, even if the contribution of order H to the mass is absent during inflation it still 
typically reappears afterwards. 

If the interaction is stronger than gravitational strength, one presumably expects a big- 
ger mass 0, but in the present context this case is not of interest because it will typically 
be associated with production mechanisms that dominate the vacuum fluctuation. In par- 
ticular, if the interaction is strong enough to maintain thermal equilibrium the mass is of 
order the temperature T, but then the thermal abundance dominates. 

We shall assume that m(t) is not varying rapidly on the timescale of the expansion rate 
H. Rapid variation will occur if m(t) is generated by a coupling to a homogeneous oscillating 
scalar field, because the oscillation frequency is necessarily big compared with H (and so 
usually is the amplitude). But in that case we are dealing with a separate mechanism of 
particle creation jpH|-^2]|, which is decoupled from the vacuum fluctuation that we consider 
here. 

With the above discussion in mind we shall assume that a species with mass m has an 
effective mass m(t) = \iH in the early universe, with fi a constant of order 1. This should 
give an order of magnitude estimate, though rj may actually be expected to have different 
values during and after inflation. 

The layout of this paper is as follows. In Section 2 we give the formalism for calculating 
the occupation number per quantum state, treating separately the cases of spin and spin 
1/2. In Section 3 we estimate the corresponding cosmological abundance and compare it 
with the abundance coming from particle collisions, on the assumption that these are of only 
gravitational strength. In Section 4 we summarize the conclusions, and some specific cases 
are worked out in an Appendix. 

II. PARTICLE CREATION 

In this section we give the basic formalism for the creation of spin zero and spin half 
particles. The reader who is interested only in the results may skip to the next section, 
where the rather simple conclusion is summarized. 
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In considering particle creation from the vacuum fluctuation, we are interested only in 
species which have very weak interactions. One reason is that a species with significant 
interactions will achieve thermal equilibrium, losing all memory of its earlier abundance. 
The other is that there are strong cosmological constraints on the abundance of a long-lived 
species, whereas short-lived species with stronger interactions decay before they can have 
any direct cosmological effect. 

In the case of a scalar field the quantum fluctuation is competing with a possible classical 
displacement of the field from its vacuum value, and this latter effect must be absent if the 
quantum fluctuation is to be significant. If the moduli fields of superstring theory exist with 
masses of order 100 GeV, a displacement for them at the classical level leads to a gross over- 
abundance in the standard cosmology, which is called the 'moduli problem' [13|-|17,6, |18 , 19 



As far as one can tell, the only way of solving the problem seems to be to dilute the moduli 
abundance with a late bout of 'thermal' inflation as discussed in detail in Refs. fT9|^2| (for 



earlier incomplete discussions of thermal inflation see the works cited in Ref. JTH|). How 



ever, this would also dilute the abundance of particles produced by the vacuum fluctuation 
that we are considering. (The vacuum fluctuation during thermal inflation is likely to be 
insignificant since the particle masses are likely to be much bigger than H, which is much 
less than 100 GeV.) Thus, for our discussion to be relevant in the context of supergravity 
the moduli problem has to be eliminated by requiring a mass ^> 100 GeV, or else no classical 
displacement, rather than by modifying the standard cosmology. 

As we are dealing with the early universe we can assume that the density is critical, so 
that the line element can be written 

ds 2 = cl 2 {ji) [drj 2 — dx 2 — dy 2 — dz 2 ] (1) 

where a is the scale factor and rj is conformal time, related to the proper time of a comoving 
observer by drj = dt/a. Using the time coordinate x° = rj and comoving space coordinates 
{x 1 , x 2 , x 3 } = {x, y, z}, the corresponding metric components are g£ = a5%. In what follows, 
a prime will denote d/dr] and an overdot will denote d/di. 

We assume that at early times there is almost exponential inflation, and set the Hubble 
parameter H = a/a equal to a constant H*. Then we can choose rj = — (aif*) -1 . After 
inflation we typically have either radiation domination when we can choose rj = (aH)" 1 or 
matter domination when we can choose rj = 2(aH)~ 1 . 



A. Spin zero particles 



The formalism in the case of a spin zero field is well known [23-25], but we reproduce it 



briefly for comparison with the less well-known spin half case. Then we apply it to the case 
m(rj) ~ H, which has not been discussed before. 

We take as the Lagrangian density of the scalar field 

C = ^(-detgy/'id^d^-im 2 ^)-^ 2 ] (2) 

We have allowed a coupling to the curvature scalar R = —6a" /a 3 = 6 [a/a + H 2 ], but will 
assume that [£[ < 1. During exponential inflation, a oc= exp(Ht) with H constant so 
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R = —12H 2 . During matter domination a oc t 2//3 so R = —3H 2 , and during radiation 
domination a oc a 1 / 2 so R = 0. 

With the exception of this coupling, we have assumed that all interactions can be repre- 
sented by a time dependent mass m(rj), and as discussed earlier we set 

m(r]) = fJ>H(ri) (fiH > m) 

m(j]) = m (/iff < m) (3) 

where m is the true mass and \x is a constant for order 1. 

In order to quantize the theory we may start with annihilation operators a^, which 
commute except for [at, a^] = 1- Then cj) is Fourier expanded in a comoving box much 
bigger than the region of interest in the form 

= a~\v) E [MvW + w* k (r])al k ] e ikx (4) 

k 

The equation of motion of w k is 

w" k + n 2 k w k = (5) 

where 

nt^e + a^w + il-OR] (6) 

o 

In the above expansion, the factor 1/a was pulled out to make Eq. (||) simple. It has the 
property that the independent solutions Wk and w k have a constant Wronskian, which we 
set equal to i, 

w k {w' k )* - w* k w' k = i (7) 
If Q k satisfies the adiabatic condition 

nyn k < n k (8) 

then w k is a sinusoidal function of rj with angular frequency Q k . With our assumptions about 
m(rj) and |£|, the adiabatic condition is satisfied whenever the expansion rate a' /a = aH 
is negligible compared with Q k . This criterion is satisfied for all all wavenumbers k in the 
late-time regime where H is less than the true mass m, whereas at earlier times it is satisfied 
only by modes with wavelength inside the horizon, corresponding to k ^> aH, 

k~>aH (H >m) 

all A; (.ff<m) (9) 

In this adiabatic regime we recover flat spacetime field theory, with Vt k /a (the physical 
angular frequency) equal to the particle energy E k = \{k/a) 2 + m 2 ] 1 / 2 . Keeping the slow 
variation of the adiabatic quantities, the canonical formulation of flat spacetime field theory 
corresponds to w k = w k at , with 
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w 



flat = 1 e -i(E k t+ Xk ) 



(10) 



where Xk is slowly varying and the slowly varying normalization follows from Eq. ([!]). These 
expressions lead to the particle concept, the number of particles per quantum state having 
the expectation value = (a^a^). 

Modes of interest are well inside the horizon now, and for inflation to do its job they 
have to be well inside the horizon at the beginning of almost-exponential inflation. (We 
assume as usual that the present density is critical, Qq = 1.) Thus they start and finish in 
the adiabatic regime, but if we choose Wk = u^ at at early times we will have at late times 
some solution Wk{jf) given by 

w k (r,) = 0^(7,) + t3 h (v$?(T,)y (11) 
From Eq. (|7|), and f3k are constants which satisfy 

K| 2 - \(3 k \ 2 = 1 (12) 
The Fourier expansion Eq. (HD can be written at late times 



w^( V )a k + w^*( V )al k \e ik - (13) 

with where = afc<2k + /^at k . Finally, we suppose that no particles exist at early times, 
(a k ak) = for all k. Then the number of particles per state at late times is (a^Ok) = |Ac| 2 . 

By solving Eq. (|5|) with the initial condition Eq. fliTf), we can read off from Eq. (|TT]). 
It is negligible, and hence there is no particle creation, if the adiabatic condition Eq. (|9]) is 
satisfied at all times. This is in general the case only for modes which are well inside the 
horizon at the end of inflation, ie., for modes with k 3> a*H* where a star denotes the end 
of inflation. For longer wavelength modes adiabaticity generally fails leading to significant 
particle creation. (Adiabaticity holds for all modes if m 2 (rj) is negligible and £ = 1/6, or 
more generally if m 2 (r]) and (| — £)R cancel, but there is no reason to expect these special 
cases to be realized.) 

In the absence of any special reason, one expects that and (3^ will be very roughly 
of the same order, which implies that very roughly they are both of order 1. One therefore 
expects that very roughly the occupation number \fik\ 2 will be of order 1 for modes which 
are inside the horizon at the end of inflation, and practically zero for shorter wavelength 
modes. In the appendix some specific cases are examined, and the expectation verified. 

It should be emphasized that we are here talking about the case where the mass is of 
order H during inflation. In the opposite case m <C H it is well known, from the case of 
axion physics that the occupation number \/3k\ 2 becomes huge in the limit of large 

wavelengths. 



B. Spin half particles 



The spin 1/2 case has been discussed briefly before especially in the approximation of 
weak particle production p7||2"8l , |i"2"! , but the following account is more systematic and explicit 
and we treat for the first time the physically interesting case m(i]) ~ H . 
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The starting point is the Dirac equation in a Robertson- Walker universe, which can be 
derived from a Lagrangian [28.12] or from more general arguments [27|. With appropriate 
normalization as specified in a moment, it is simply the ordinary Dirac equation written in 
comoving coordinates, 

(za~ -m + £hj°)iI) = (14) 

To quantize one can start with the usual anticommuting annihilation operators a rk and 6 rk 
(where r denotes the spin state), and make the Fourier expansion 

$ = a" 3 / 2 £ £ [ Mrk a rk + v rk bl_^] e- ik ' x (15) 

k r=±l 

(This applies to a Dirac field; for a Majorana field a rk = 6 rk which will make no difference 
in what follows.) Substituting this expression into the Dirac equation and choosing for each 
mode the z axis along k gives 

ijou' rU - [j z k + am]w rk = (16) 
hoK-k - hzk + am]f r _ k = (17) 

Any solution of this equation can be written as a linear combination of Wi k , «2k, v i-k an d 
We impose on these four objects the usual orthonormality condition 

M Ik^k = vl k v sk = 5 rs (18) 
uJk^k = (19) 

This condition is preserved by the evolution equation Eq. (14), because (7 ) 2 is the unit 
matrix and 7°(a _1 7.k + m) is Hermitian. Note that we did not pull out a factor l/(2Ej s .) 
before defining u k and t> k , as is often done in flat spacetime. Had we done so, w k and v k 
would have satisfied a more complicated evolution equation because of the time dependence 
of Ek- We did however pull out a factor a -3 / 2 , which again is crucial in making the equation 
simple. 

As the evolution of the two spin states is independent we need consider only one of them. 
This means that we can regard if) as only a two component object so we drop the subscript 
r. The Dirac matrices can be taken to be 

-(:_ o *-{y ) <*> 

Now drop the subscripts k and — k, and denote the upper and lower components of u 
and v by a subscript ±. They satisfy uncoupled equations 

u'± + a 2 [El ± iHm]u± = 

v'l + a 2 {E 2 k ± iHm]v± = (21) 
where as before E^ is the energy related to the momentum = k/a by E\ = m(r]) 2 + p 2 . 
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With our assumption about m(j]), Hm <C E\ during any era in which the adiabatic 
condition Eq. (J|) is satisfied, and we recover flat spacetime field theory with the canonical 
choice 

Jat_ ( E k + m \^( 1 \ i{Ekt+Xk) 

^flat = ( ^+ W ) 1/2 ^ -Pfc/(^ fc + m) ^ ei{Ekt+Xk) (22) 

The ratio of the upper and lower components comes from Eq. ([H]), Xk is slowly varying, 
and the slowly varying normalization comes from conditions ( [L8D and flT9l). As one knows 
from flat spacetime field theory, these expressions lead to the particle concept, the number 
of particles per quantum state being (a^a^) and the number of antiparticles being (&k^k)- 
Just as for a scalar field, we consider solutions u and v, which become equal to w flat and 
^flat we jj b e f ore horizon exit during inflation. Remembering that u + has the same evolution 
equation as u*_ and similarly for v±, they are of the form 



u = 

\Ji ) 

(23) 

Thus, in the late time adiabatic regime they are of the form 

u(r ] ) = a k u &at (r ] )+f3 k v Rat (r ] ) (24) 
v{ V ) = alv^\ V )-Plu^\ V ) (25) 




Here a k and (3 k are constants, which from the normalization condition ( |T8| ) satisfy 

\a k \ 2 + \(3 k \ 2 = 1 (26) 
The contribution of the mode to ip is 

a" 3 / 2 (u Rat ~a k + t; flat 6l k ) (27) 

where a k = a k a k — /3^&! k and 6 k = /5fcak + a^fel k . Finally, we suppose that no particles exist 
at the beginning of inflation, (a^ak) = (b^b^) = 0. Then the number of particles per state 
at late times is (fik^k) = \f3 k \ 2 , and the number of antiparticles (&k^k) is the same. 

In order to calculate \f3 k \ 2 precisely one has to solve the second order differential equation 
(^T|) with the initial condition (|l0t). But since it must lie between and 1 one may hope 
that an adequate order of magnitude estimate will be obtained by setting it equal to 1 if 
there is an era of non-adiabaticity and equal to otherwise. A couple of cases are worked 
out in the Appendix, and they support this hope provided that /i is somewhat less than 1. 
(For fi ~ 1 one finds \f3 k \ 2 ~ e~ 27TtI .) 

As in the spin zero case one will have adiabaticity at all times for those modes which are 
far inside the horizon at the end of inflation. In addition, there is adiabaticity at all times 
if the mass is negligible in the sense that p k (r]) ^> m(i]) and p(rj) 2 3> Hm(r]), because then 
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Eq. (pip reduces to u± + k 2 u± = 0. However if, as we are supposing, 777,(77) ~ ^ (until the 
epoch H ~ m, which we are supposing occurs after the end of inflation), these conditions 
are equivalent to the adiabatic condition Pk{f]) ^> H. The conclusion is therefore that if 
771(77) ~ H (at least during inflation) one may hope for significant particle production in all 
modes for which the adiabatic condition fails just as in the spin case. 

III. COSMOLOGICAL IMPLICATIONS 

We have considered a particle, whose true mass m is less than the Hubble parameter 
during inflation. We have assumed that its mass 771(77) m the early universe, before the 
epoch H ~ m, is of order H and smoothly varying which is what one expects in general 
for a particle with an interaction of only gravitational strength. On these assumptions, we 
have found that of order 1 particle per quantum state is created for modes which are outside 
the horizon at the end of inflation (corresponding to k < a*H*), whereas for smaller scales 
particle creation is insignificant. Although we have demonstrated this only for spin zero 
and spin half particles the argument looks quite general and it reasonable so suppose that 
it holds for any spin.[] 

In this section we see what the estimate implies for cosmology. We assume of course that 
the species does not thermalize subsequently, since then it loses all memory of its earlier 
abundance. 

A convenient measure of the abundance is provided by the ratio n/s, where n is the num- 
ber density and s is the entropy density. After inflation is over and reheating has occurred, 
this ratio can be taken to be constant between the epoch when the particle stops interact- 
ing and the epoch when it decays. For a species which is initially in thermal equilibrium, 
n/s ~ g~ 1 where g* ~ 10 to 10 3 is the effective number of particle species at the epoch when 
thermal equilibrium fails. 

Summing over all modes, the vacuum fluctuation gives 

77 = (4tt 2 )-V 3 \(3 k \ 2 k 2 dk ~ \(3 a * H »\ 2 {a*Hjaf (28) 

Jo 

In estimating the integral we have assumed that it is dominated by the upper limit. This is 
true provided that \Pk\ 2 ~ k 1 with 7 > —3, which holds for the cases we evaluated in the 
last section. As discussed in the text we expect \f3 a *H*\ ~ 1 and from now on we adopt the 
value 1. 

The energy density just after reheating is of order gT% where Tr is the reheat tem- 
perature, and the entropy density is of order gT^. Suppose first that reheating is prompt 
so that gTft ~ V* where V* = 3iJ^M|, ; is the potential at the end of inflation. (Here 
Mpi = (8ttG)" 1/2 = 2 x 10 18 GeV is the reduced Planck mass.) Then 



As mentioned earlier, the gravitino is practically massless during inflation according to some 
models but it will have the canonical mass ~ H afterwards which should be enough to ensure that 
the stated result holds. If the mass were much less than H even after inflation one would need 
to ask whether the gravitino field equation is conformally invariant since in that case gravitino 
production would be negligible. We are not sure of the answer to that question. 
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n/s ~ {HjM Pl f 2 ~ {V^/M Pl f (29) 
If reheating is delayed this estimate is reduced by a factor Tr/U 1 / 4 to become 



n / s ~ ~T~r~ ( "T7 — I (30) 
7 M Pl [M Pl ) 1 ; 

Now let us compare this abundance with that coming from particle collisions, assuming 
that the interaction of the species under consideration with other particles is of only grav- 
itational strength. Candidates for this case are the gravitino (spin 3/2) the moduli of 
superstring theory (spin 0) and their superpartners the modulini (spin 1/2) flI5,IE|l. Up to 



numerical factors and dimensionless couplings, the cross section for producing the species in 
a typical relativistic particle collision process is of order Mp ( 2 , and as a result the abundance 



arising from this mechanism is [30 



n/s = JTr/Mpi (31) 

where / comes from the numerical factors and couplings. In the case of a gravitino (with 
the standard properties), / can be estimated in a fairly model- independent fashion [|(J to 
be of order 10~ 5 . It is not quite clear what value of / to expect more generally, for example 
for moduli and modulini. 

Comparing these two estimates, we see that creation from particle collisions dominates 
creation from the vacuum fluctuation, provided that V 1 ^/M P i < f 1 ^ 2 . In order not to create 
too much cmb anisotropy one must have |I| yl,i < 10~ 2 M Pi , and we see that if the upper 



bound is attained the vacuum fluctuation mechanism is marginally stronger, at least for 
gravitinos. 

Now let us consider the cosmological constraints. First suppose that the species is un- 
stable. Then it should decay with only gravitational strength or as we have just discussed 
the abundance from the vacuum fluctuation cannot be significant. The decay rate is then 
T = 10~ 2 7m 3 /Mp Z where the coefficient 7 is expected |L9] , ^lJ to be < 1. The correspond- 
ing decay time T _1 is typically after the beginning of nucleosynthesis, and as a result the 
abundance during nucleosynthesis is constrained to be |30| 



n/s < 1(T 12 to HT 15 (32) 

Taking n/s < 10~ 15 , the constraint on the reheat temperature following from the abundance 
due to particle collisions is 

Tr < 10- 15 f- l M Pl (33) 

Assuming that (standard) gravitinos exist they correspond to / ~ 10 -5 which implies the 
well known constraint Q 

T R £l0 8 GeV (34) 

Using the vacuum fluctuation abundance instead multiplies this result by a factor 
/(Mp;/^ 1 / 4 ) 2 , which as we just noted might by marginally bigger than 1. 
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Now consider a stable particle species. The requirement Oo < 1, where Q,q is the present 
density in units of critical density, gives |l[ 

n/s<(leV/m) (35) 

where m is the mass of the species. Using the abundance from particle collisions, this gives 

/ m \ /10 5 \ Zl0 16 GeV\ 

(lOMiv) S (— )HH (36) 
Using instead the abundance Eq. ( |30D coming from the vacuum fluctuation it gives 




10 5 GeW ~ V yl/4 
If reheating is instantaneous this becomes 

10 16 GeV\ 3 



1 MeV/ ~ 1 V 1 '* ) 



IV. CONCLUSION 

We have presented a complete and unified formalism for the cases spin zero and spin 
half. It invokes minimal physical assumptions, so that for example canonical quantization 
is required only in the adiabatic regime, with the classical evolution of the (Heisenberg 
representation) operators being the only quantum physics that is used in the intervening 
regime. The formalism is more complete than earlier discussions for the spin half case, and 
in both cases it is here applied for the first time to the case of an effective mass ~ H . 

The bottom line consists of the bounds on the parameters presented in the previous 
section, associated with the existence of particles having gravitational strength interactions. 
Though they are are at best only comparable with similar ones derived on the assumption 
that the particles are created through collisions after reheating, they are of interest in prin- 
ciple, and might be relevant if the inflation scale is high and production through collisions 
is more suppressed than for the gravitino. 

APPENDIX 

In this appendix some specific cases are worked out. 

Spin zero 

For modes which are well inside the horizon at the epoch H = m (equivalently, for modes 
which are relativistic at that epoch), it is enough to solve Eq. (|5[) with m(rj) = fiH. It then 
has the form 
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k 



w k = 



(39) 



During exponential inflation, radiation domination and matter domination, this is equivalent 
to a Bessel equation for Wk/f] 1 ^ 2 , with the order v given in Table 1, and assuming that the 
universe makes abrupt transitions between these regimes we can calculate (3k by matching 
Wk and its first derivative at each transition. 

Let us drop the subscript k, and let u> in f denote w during inflation. For argument x 1, 
the Bessel function becomes 



(40) 



By virtue of our assumptions, tf nf = | — fi 2 — 12£ is not much bigger than 1 so the regime 
of large argument corresponds precisely to the early-time adiabatic regime krj ^> 1. In this 
same regime Ef. = k/a, so one sees that the solution with the required early-time behaviour 
Eq. (|TDp is 



Winf 



^(-V^H^i-kv) ex P [(2^ + l)^z] 



(41) 



We shall assume that inflation gives way quickly to an era of either radiation or matter 
domination, which persists until the adiabatic regime is encountered. Then, the solution 
w evSu of Eq. (|5|) which reduces to Eq. ([H]) at late times is the solution w rac \ which has the 
required form Eq. (|i~0|) at late times is 



Wo 



71 ' ^H^i-kr,) exp[-(2z/ cra + 



(42) 



2 ' "era \ '/ * L \ / ^ 

We are looking for a solution of the form 

w = aw eia + (3w* er& (43) 

whose value and first derivative matches the inflationary solution. Remembering the nor- 
malization Eq. (|7|) one has one finds 

a = i(w ini w* e ^ - w' in{ w* CVSL ) 

(3 = i(w inf w' era - W^era) (44) 

For modes with k 3> a^H* the evolution is adiabatic and \(3\ is negligible. Considering 
for simplicity only the opposite case k <C a*H*, we can use the small small-argument limits 
of the Bessel functions, 

-ill 



H?-\x) 

Hl 2 \x) 
Using them one finds 



sini/vr \T(1 - v) \2j T(l + u) \2 
+i ( 1 fx 



sin z/7r \ r( 1 — v) V2 



T(l + u) V2 



(45) 
(46) 



P = — (X(z/ inf , ^ era ) + X(-V iv£ , ~Z/ era ) + X(z/ inf , -V eTa ) + X(-Z/ inf , 2/ era )) (47) 
47T 



where 



X(a, 6) = (1 + a + 6)e^" a )f r(-a)r(-6) (?^\ 



a+b 



(48) 
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Radiation domination after inflation 

For the case of a scalar field with conformal coupling (£ = |) the equations of motion 
during radiation domination are identical to those during inflation (ie z/ in f = f ra d)- For this 
situation the equation for (5 simplifies significantly 

/^ = if(i + 2-)r(-,) 2 fi-4rl" + (i-2-)r(,) 2 fi- 1 S "' 



4:7i \ \2a*H*J \2a*H* 

2tt cos 7TlA 



v Sm7Tf 



(49) 



For the case where c = /iwe have v = i^J~^ and hence the third term dominates 

% cosh \ %TX 

(3 ~ (50) 

y/3 sinh J jTT 

\P\ 2 *l (51) 

In general the other two (neglected) terms in the formula for f3 will produce oscillatory terms 
in the particle numbers, but do not affect the total number significantly. 

For the case of a scalar field without coupling to the Ricci scalar v will in general be 
different during the different era's under consideration. During inflation v = \J\ — (J> 2 and 

during radiation domination v = \J\ — /i 2 - For \ < /i < | then z/ in f will be real whilst z/ ra d 
will be imaginary. In this case one of the terms in the expansion will dominate. If we set 
z/ in f = a and z/ rad = ib for convenience, where both a and b are real, then taking the leading 
term in our expression gives 

(1 h \ ~ a ~* 
(52) 

\[3\ 2 ~ , - - ((1 - a) 2 + 6 2 )r(a) 2 (53) 

IPI 167r6sinh7r6 VV ; J K J \2a*Hj K ' 

In the case where n = 1 this works out as 

w^HItjl)^ (54) 



Matter domination after inflation 

If a conformally coupled scalar field re-enters the particle horizon during during matter 
domination as opposed to radiation domination as discussed above then the equations de- 
scribing its evolution will no longer be the same during the two epochs and hence the full 
general solution must be used. Taking a particular choice of mass such that fi = 1 then the 
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orders of the two bessel functions during the subsequent epochs will be u in f = i a/3/4 during 



inflation and v mat = ly 15/4 during matter domination. For these values it is clear that the 
third term in the general expression will dominate, and hence our approximate solution is 

/ i \ -il-07 

1 / - - \ 1 A„„ I . /„ , . \ „ / . /_ _ , . \ / rC 



^ -(1 - 1.07,)e-r ^3/4j r ^ j (55) 

m 2 - ^^ 8 1 r (^) H r (iy/W) I 2 (56) 

« 0.3 (57) 

For a minimally coupled field, £ = 0, the equation for /3 is the same as that for re- 
entry during radiation domination with the exception that the constant b is replaced by the 
corresponding value for matter domination (ie v mat = ib) . Calculating the result for the case 
where fi = 1 gives 

w - 003 y^r <58 » 



Spin half 

In contrast with the spin zero case, we definitely need to consider only modes which are 
well inside the horizon at the epoch H = m, corresponding to k ^> a*H*. The reason is 
that modes with smaller wavenumber are killed by the phase space factor k 2 in Eq. (pgp, 
owing to the fact that the occupation number can can never exceed 1. Thus it is enough 
to set m(r]) = fiH at all times. Then, during exponential inflation, matter domination and 
radiation domination Eq. ([H|) is equivalent to a Bessels equation for u±/i] 1 ^ 2 and v±/ r q 1 ^ 2 
with the index given in Table 1. During inflation, the solutions reducing to Eq. ( f22D at early 
times are 

/inf(i) = l^^\e^ H^(k\ V \) (59) 

/inf(2) = -'-J^le-^H^JkW) (60) 

We have used this formalism to estimate the occupation numbers on the assumption 
that inflation is promptly followed by either radiation or matter domination. The solutions 
during radiation domination which reduce to Eq. (^) at late times are 

/radd) = -~ ^e-^Hf_ ifi (k V ) (61) 

/rad(2) = \ y/^kqe^HPjkrj) (62) 

The solutions for matter domination are the same with /i replaced by 2/i. Matching either 
set of solutions gives 
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And hence 



/inf(l) = a/era(l) _ Pfeia,(2) 
/inf(2) = a /cra(2) + /^/era(l) 



P — (/kif(2)/era(l) ~ /inf(l)/era(2)) 



where the right hand side is evaluated at the end of inflation. 
For radiation domination this gives 



P 



1 

2^ 



i it 



2ifi 



2 ciji. 



— c.c. 



and for matter domination it gives 



(63) 
(64) 



(65) 



(66) 



where 



13 = — (Be^ - B*e~^) 



B = r(i + W r ( I + 2M.)(-^ 



-3ifj, 



(67) 



(68) 



Since |T(| ±iy)| 2 = 7r/cosh(7ry) these expressions give \(3\ 2 < 1 as required. In the 
adiabatic regimes k 3> a*if* and /i 3> 1 they give the expected result |/3| 2 <C 1. In the 
regime k <C a^H* and fi ~ 1, they make |/3| 2 a rapidly oscillating function of k with a mean 
value roughly of order e~ 2nfl , which is roughly of order 1 with our assumed value /i ~ 1. 



TABLE 1. Tl 


le order of the Bess< 
Scalars 


A function 
Spinors 


Inflation 
Radiation dom. 
Matter dom. 


z, 2 = § - 12£ - 4/i 2 


v — \ — i/j, 
v — \ — i/i 
v — | — 2i\x 



14 



Acknowledgements 

We are indebted to Ewan Stewart for many helpful discussions about supergravity. The 
work is partially supported by grants from PPARC and from the European Commission 
under the Human Capital and Mobility programme, contract No. CHRX-CT94-0423. 



15 



REFERENCES 



[1] E. W. Kolb and M. S. Turner, The Early Universe, Addison- Wesley, New York (1990). 
[2] A. D. Linde, Particle Physics and Inflationary Cosmology, Harwood Academic, Switzer- 
land (1990). 

[3] M. Dine, W. Fischler and D. Nemeschansky, Phys. Lett. 136B, 169 (1984). 
[4] O. Bertolami and G. G. Ross, Phys. Lett. B183, 163 (1987). 

[5] E. J. Copeland, A. R. Liddle, D. H. Lyth, E. D. Stewart and D. Wands, Phys. Rev. 
D49, 6410 (1994). 

[6] M. Dine, L. Randall and S. Thomas, Phys. Rev. Lett. 75, 398 (1995). 
[7] D. H. Lyth, |hep-ph/9609431 . 



E. D. Stewart, Phys. Rev. D51, 6847 (1995); E. D. Stewart, Phys. Lett. B345, 414 
(1995). 

[9] E. D. Stewart, personal communication. 
[10] J. H. Traschen and R. H. Brandenberger, Phys. Rev. D42, 2491 (1990); L. Kofman, 
A. D. Linde and A. A. Starobinsky, Phys. Rev. Lett. 73, 3195 (1994); Y. Shtanov, J. 
Traschen and R. Brandenberger, Phys. Rev. D51, 5438 (1995); D. Boyanovsky et al., 
Phys. Rev. D51, 4419 (1995); |hep-ph/9505220| ; |hep-ph/95074T4] ; |hep-ph/951136T| ; [Kep] 



ph/9602232] ; |iep-ph/9603"2T7| ; |hep-ph/9606205| ; |hep-ph/9608"208| ; I. I. Tkachev, Phys. 



Lett. B376, 35 (1996); M. Yoshimura, |hep-th/9506T76| ; D. I. Kaiser, |astro-ph/950TT08| ; 
L. Kofman, A. D. Linde and A. A. Starobinsky, Phys. Rev. Lett. 76, 101 (1996); D. 
T. Son, |hep-ph/ 960 1377| ; A. Riotto and I. I. Tkachev, |hep-ph/9604l4l] ; E. W. Kolb, 
A. Linde and A. Riotto, |hep-ph / 9606260j ; G. W. Anderson, A. Linde and A. Riotto, 
|hep-ph/96054T^ ; G. Dvali and A. Riotto, |hep-ph/ 960M3"Tt A. Riotto, E. Roulet and I. 
Vilja. [hep-ph/9607403| . L. Kofman, A. Linde and A. Starobinsky, |hep-ph/ 9608341 . 

[11] A. D. Dolgov and D. P. Kirilova, Soc. J. Nucl. Phys. 51, 273 (1994). 

[12] A. D. Dolgov and K. Freese, Phys. Rev. D51, 2693 (1995). 

[13] G. D. Coughlan et al., Phys. Lett. 131B, 59 (1983). 

[14] J. Ellis, D. V. Nanopoulos and M. Quiros, Phys. Lett. B174, 176 (1986); B. de Carlos, 

J. A. Casas, F. Quevedo and E. Roulet, Phys. Lett. B318, 447 (1993). 
[15] T. Banks, D. B. Kaplan and A. E. Nelson, Phys. Rev. D49, 779 (1994). 
[16] T. Banks, M. Berkooz and P. J. Steinhardt, Phys. Rev. D52, 705 (1995). 
[17] L. Randall and S. Thomas, Nucl. Phys. B449, 229 (1995). 

[18] P. Binetruy and M. K. Gaillard, Phys. Rev. D 34, 3069 (1986); G. Lazarides and Q. 
Shafi, Nuc. Phys. B392, 61 (1993); D. H. Lyth and E. D. Stewart, Phys. Rev. Lett. 75, 
201 (1995). 

[19] D. H. Lyth and E. D. Stewart, Phys. Rev. D53, 1784 (1996). 
[20] T. Barriero et al., |hep-ph/9602263| . 

[21] E. D. Stewart, M. Kawasaki, and T. Yanagida, |hep-ph / 9603321 . 

[22] K. Choi, E. J. Chun and E. J. Kim, |hep-ph/9608222| . 

[23] For a review, see A. R. Liddle and D. H. Lyth, Phys. Rep. 231, 1 (1993). 

[24] T. Damour and A. Vilenkin, |hep-th/950314~9| . 

[25] T. Damour and A. Vilenkin, Phys. Rev. D53,2981 (1996). 

[26] D. H. Lyth, Phys Rev D 45 3394-3404 (1992). 

[27] L. Parker, Phys. Rev. D3, 346 (1971). 

[28] S. G. Mamaev, V. M. Mostepanenko and V. M. Frolov, Sov. J. Nucl. Phys. 26, 1 



16 



(1977); N. D. Birrel and P. C. W. Davies, Quantum Field Theory in Curved Space- 
time, Cambridge University Press (1982). 
[29] For reviews of supersymmetry, see H. P. Nilles, Phys. Rep. 110, 1 (1984) and D. Bailin 
and A. Love, Supersymmetric Gauge Field Theory and String Theory, IOP, Bristol 
(1994). 

[30] J. Ellis, J. E. Kim and D. V. Nanopoulos, Phys. Lett. 145B, 181 (1984); J. Ellis et al, 

Nucl. Phys. B373, 399 (1992). 
[31] D. H. Lyth, Phys. Lett. 147B, 403 (1984); D. H. Lyth, Phys. Lett. 150B, 465 (1985); 

D. H. Lyth, Phys. Rev. D 31 (1985) 1792. 



17 



